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Abstract 

A  study  is  conducted  on  a  cogeneration  system  that  incorporates  a  natural  gas  fed  internal-reforming  solid  oxide  fuel  cell  (IRSOFC)  and  a 
zeolite/water  adsorption  chiller  (AC).  The  main  aim  is  to  investigate  the  performance  of  this  combined  system  under  different  operating  conditions 
and  design  parameters.  A  mathematical  model  is  developed  to  simulate  the  combined  system  under  steady- state  conditions.  The  effects  of  fuel 
flow  rate,  fuel  utilization  factor,  circulation  ratio,  mass  of  adsorbent  and  inlet  air  temperature  on  the  performance  are  considered.  The  results  show 
that  the  proposed  IRSOFC-AC  cogeneration  system  can  achieve  a  total  efficiency  (combined  electrical  power  and  cooling  power)  of  more  than 
77%.  The  electrical  efficiency  is  found  to  decrease  as  the  fuel  flow  rate  increases,  while  the  cooling  power  increases  to  a  constant  value.  The  total 
efficiency  reaches  a  maximum  value  with  variation  of  the  fuel  utilization  factor.  Both  the  circulation  ratio  and  the  inlet  air  temperature  exert  positive 
impacts  on  system  efficiency. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  have  been  considered  to  be 
promising  energy  conversion  technologies  because  they  can 
provide  a  highly  efficient  rate  of  energy  conversion  with  low 
pollutant  emissions.  The  chemical  energy  of  the  fuel  gas  is  con¬ 
verted  directly  to  electricity  in  a  SOFC,  and  hence  it  is  expected 
that  high  electrical  efficiencies  can  be  achieved.  The  high  oper¬ 
ating  temperature  of  SOFCs  allows  the  use  of  cogeneration 
and  hybrid  systems.  Gas-solid  adsorption  cooling  systems  are 
environmental-friendly  compared  with  traditional  CFC  systems 
based  on  the  vapour  compression  cycle  as  they  employ  safe 
and  non-polluting  refrigerants.  Another  advantage  of  adsorp¬ 
tion  cooling  systems  is  that  they  can  be  driven  by  waste  heat 
or  solar  energy.  As  a  result,  both  SOFCs  and  adsorption  cool¬ 
ing  systems  have  attracted  much  research  attention  in  the  recent 
years. 

Fuel  cell  cogeneration  systems  (FCCS)  have  been  applied 
successfully  in  some  developed  countries.  There  have  been 
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increased  interests  in  such  systems  because  they  can  achieve 
higher  efficiency  in  electricity  production  and  a  higher  overall 
efficiency  with  a  low  level  of  adverse  environmental  impact. 
One  type  of  fuel  cell  and  absorption  cogeneration  system  has 
been  investigated  by  several  researchers  [1,2]  to  produce  elec¬ 
trical  power  and  cooling  simultaneously.  Because  the  gas-solid 
adsorption  cooling  system  is  more  flexible  in  operation  com¬ 
pared  with  the  absorption  cooling  cycle  [3],  a  combined  fuel  cell 
and  adsorption  cooling  cogeneration  system  can  be  a  promising 
alternative. 

In  this  paper,  an  internal-reforming  solid  oxide  fuel  cell 
(IRSOC)  and  a  zeolite/water  adsorption  chiller  (AC)  cogenera¬ 
tion  system  is  proposed.  A  numerical  model  based  on  the  first 
law  of  thermodynamics  is  developed  to  simulate  this  combined 
system  under  steady-state  operating  conditions.  The  main  aim 
is  to  investigate  the  performance  of  the  combined  system  under 
different  operating  conditions  and  parameters. 

2.  System  description 

A  schematic  of  the  proposed  SOFC  and  adsorption  chiller 
cogeneration  system  is  shown  in  Fig.  1 .  The  system  incorporates 
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Nomenclature 

a  flow  rate  (kmol  h- 1 ) 

E  cell  voltage  (V) 

Evq  reversible  cell  voltage  (V) 

F  Faraday  constant  (96,487  C  mol- 1 ) 
i  current  density  (mA  cm-2) 

io  exchange  current  density  (mA  cm-2) 

I  current  (A) 

Kv  equilibrium  constant 

n  number  of  electrons  participating  in  the  reactions 

p  pressure  (bar) 

Q  heat  transfer  rate  (kJ  h- 1 ) 

ra  circulation  ratio 

R  universal  gas  constant  (8.314  J mol-1  K-1) 

T  temperature  (K) 

Uf  fuel  utilization  factor 

Greek  letters 

/3  transfer  coefficient 

8  thickness  (cm) 

AH  enthalpy  of  formation  (kJ  mol-1) 

rjact  activation  over-potential  (V) 

Oconc  concentration  over-potential  (V) 

ohmic  ohmic  over-potential  (V) 

rjt  total  efficiency 

Subscripts 
a  anode 

c  cathode 

r  reforming  reaction 

s  shifting  reaction 

Superscripts 
f  fuel 

i  inlet 

o  outlet 


a  natural  gas  feed,  an  IRSOFC,  an  AC,  a  dc-ac  inverter,  and 
two  heat  exchangers.  The  system  can  produce  electricity,  cool¬ 
ing  and  heating  simultaneously.  Air  is  preheated  at  the  heat 
exchanger  before  it  enters  the  IRSOFC  generator.  Part  of  the 
exhaust  fuel  gas  is  recirculated  and  mixed  with  fresh  fuel  gas 
and  the  mixture  re-enters  the  IRSOFC.  Fuel  gas  and  air  flow 
through  the  anode  and  cathode,  respectively.  Fuel  gas  is  inter¬ 
nally  reformed  at  the  anode  and  hydrogen  is  produced.  The 
electrochemical  reaction  occurs  in  the  SOFC  and  electricity  is 
generated  together  with  heat.  Subsequently,  the  non-recirculated 
portion  of  the  effluent  fuel  gas  and  depleted  air  flow  into  a  com¬ 
bustor.  After  reacting  in  the  combustor,  the  fuel-air  mixture 
enters  heat  exchanger  1  where  the  air  is  preheated.  The  com¬ 
bustion  gases  exiting  from  heat  exchanger  1  then  flows  through 
the  adsorption  chiller  to  drive  this  cooling  system  during  the 
heating  process.  Zeolite/water  is  used  as  the  working  pair  for 
the  adsorption  chiller.  Zeolite  has  satisfactory  adsorption  ability 


Fig.  1.  Schematic  of  SOFC  and  adsorption  chiller  cogeneration  system. 


and  its  adsorption  properties  are  also  stable  under  high  temper¬ 
ature.  Finally,  the  gases  from  the  adsorption  chiller  are  passed 
through  heat  exchanger  2  for  heat  generation  and  then  exhausted 
to  the  environment. 


3.  Numerical  modelling 

The  numerical  model  consists  of  three  parts:  (i)  inter¬ 
nal  reforming  and  electrochemical  reaction  model;  (ii)  SOFC 
model;  (iii)  adsorption  cooling  cycle  model.  The  following 
assumptions  are  made:  (i)  equilibrium  reforming  and  shifting 
reactions;  (ii)  the  exit  temperature  of  the  cathode  is  equal  to  that 
of  the  anode;  (iii)  the  SOFC  and  heat  exchangers  are  in  thermal 
balance;  (iv)  the  burner  efficiency  is  assumed  to  be  100%;  (v)  the 
mass  of  the  adsorbent  is  proportional  to  the  surface  area  of  the 
adsorber;  (vi)  there  is  no  heat  loss  in  the  cogeneration  system. 


3.1.  Internal  reforming  and  electrochemical  reaction  model 


The  recirculated  anode  exhaust  fuel  gas  and  the  high- 
temperature  conditions  greatly  assist  the  internal  reforming  and 
shifting  reactions  in  a  SOFC.  The  mechanisms  for  internal 
reforming  and  electrochemical  reactions  for  the  SOFC  are  com¬ 
monly  reported  in  the  literature  [4-7].  The  reactions  are: 
reforming  reaction: 

CH4  +  H20  ->  CO  +  3H2  (1) 

shifting  reaction: 

CO  +  H20  ->  C02  +  H2  (2) 

electrochemical  reaction: 


H2  +  l/202  — >  H20 


where  the  reforming  and  shifting  reactions  are  assumed  to 
reach  thermodynamic  equilibrium  with  the  respective  equilib¬ 
rium  constants  represented  by: 


Pr2  •  Pco 
Pcu4  •  Pu2o  ’ 


P co2  •  P h2 
Pco  •  Pu2o 
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The  equilibrium  constants  are  considered  to  be  functions  of 
temperature  and  can  be  obtained  from  the  following  equation: 

log  Kp  =  AT4  +  BT3 +CT2  +  DT  +  E  (5) 


The  values  of  the  parameters  in  Eq.  (5)  have  been  reported  in 
literature  [8].  It  is  assumed  that  the  fresh  fuel  gas  is  composed  of 
CH4,  N2  and  CO2.  The  equilibrium  constants  can  be  rewritten 
as: 


f 

\  a\+2x  ) 


(aco,+y\  2 

y  a\+2x  J  ^cell 


(  x~y  ^  ( z~x~y^\ 
\a\+2x  )  \  a\+2x  ) 


where  x,  y  and  z  are  the  reactant  molar  flow  rates  of  CH4,  CO2 
and  H2,  respectively,  and  the  superscript  f  denotes  fresh  fuel  gas. 

The  reactant  flow  rate  z  of  hydrogen  can  be  obtained  from 
the  fuel  utilization  factor  as  follows: 


Z  —  +  aU2  +  aCo) 


where  the  superscript  i  represents  inlet  cell  and  Uf  is  the  fuel 
utilization  factor. 

If  the  value  of  the  circulation  ratio  ra  is  known,  Eq.  (8)  can 
be  rewritten  as: 


4Uf  •  afcll4 
1  -  ra  +  ra  •  Uf 


3.2.  SOFC  model 


The  cell  voltage  calculation  is  the  core  of  the  SOFC  model. 
The  cell  electrical  power  can  be  calculated  as  the  product  of 
the  cell  current  and  the  cell  voltage,  while  the  cell  current  can 
be  easily  obtained  from  a  known  value  of  the  fuel  utilization 
factor.  Based  on  the  work  of  several  researchers  [6,9,10],  the 
cell  voltage  is  given  by: 

E  =  EXq  —  (?7act  T  Oohm  T  Oconc )  (10) 


where  Ere  is  the  ideal  reversible  cell  voltage  and  the  other  three 
terms  represent  the  activation,  ohmic  and  concentration  overpo¬ 
tentials. 

The  ideal  reversible  cell  voltage  can  be  evaluated  from  the 
Nemst  equation,  i.e. 


RT 

2F 


In 


(  ^h2Q 
\Ph2P<£) 


RT\nKp,Q 
2  F 


RT  (  Ph2o 
In  1 


IF 


\ 


1/2 

P^Po2 


(11) 


where  KV  Q  is  the  equilibrium  constant  of  the  electrochemical 
reaction. 


Adsorption 

unit 


Fig.  2.  Schematic  of  adsorber  in  adsorption  chiller. 


The  anode  and  cathode  activation  overpotentials  can  be 
described  by  the  Bulter-Volmer  equation  [11]: 


exp  (1 


(12) 


where  /3  is  the  transfer  coefficient  and  iq  is  the  apparent 
exchange-current  density. 

The  ohmic  overpotential  as  suggested  by  Bessette  [11]  is 
given  by 

^ohmic  —  —  i'}  "jAj  exP  ^ j  (13) 

j  j  ^  ' 


The  parameters  Aj  and  Bj  are  constants  for  different  compo¬ 
nents,  while  8j  is  the  thickness  of  the  cell  components. 

The  concentration  overpotential  is  a  result  of  the  reac¬ 
tant  diffusion  during  reaction.  According  to  Bove  et  al.  [10], 
concentration  losses  can  be  ignored  under  normal  operating 
conditions. The  SOFC  current  can  be  quantified  using  the  fuel 
utilization  factor,  Uf,  i.e., 

I  =  ZnF  =  2FUf{4alQY{A  +  ^h2  +  aco )  (I4) 

The  SOFC  only  converts  part  of  the  chemical  energy  of  the 
fuel  into  electrical  power  while  the  remainder  will  become  heat 
to  increase  the  temperature  of  the  outlet  effluent  gas.  From  the 
following  energy-balance  equation,  the  temperature  of  outlet  gas 
is  computed. 


y>‘cA//‘c+y>°aA//i 


a 


J 


=  EAA^c  +  E<aAH°a  +  £/ 


(15) 


As  the  temperature  is  known  from  Eq.  (15)  and  for  a  specified 
fuel  utilization  factor  and  inlet  conditions,  the  species  concen¬ 
tration  in  the  SOFC  can  be  obtained  by  solving  Eqs.  (6)-(8). 


3.3.  Adsorption  cooling  cycle  model 

After  preheating  the  air,  the  gas  from  combustor  is  directed 
to  the  zeolite/water  adsorption  chiller.  The  adsorber  consists  of 
a  number  of  adsorption  units,  as  shown  schematically  in  Fig.  2. 
Every  adsorption  unit  is  a  metal  tube  covered  with  a  zeolite 
adsorbent  bed.  The  configuration  is  the  same  as  that  of  the  adsor¬ 
ber  described  in  our  previous  paper  [12].  The  exhaust  gas  from 
heat  exchanger  1  flows  through  the  metal  tube  of  the  adsorption 
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unit  and  heats  up  the  zeolite  adsorbent  bed.  The  heat  and  mass 
balance  equations  for  every  adsorption  unit  here  are  the  same 
as  those  reported  previously  [12]  except  that  the  heat-exchange 
fluid  is  changed  to  the  exhaust  gas.  Adsorption  chillers  with  a 
continuous  adsorption  cycle  are  employed  in  this  system.  Hence, 
the  cooling  power  can  be  presented  as: 


N  •  msAq[L(Te)  -  Cpl(rc  -  Te)] 

hi 


(16) 


where  N  is  the  number  of  adsorption  units;  ms  is  the  mass  of 
the  adsorbent  for  every  adsorption  unit;  th  is  the  duration  of  the 
heating  process;  Aq  is  the  change  of  adsorbed  water  uptake. 


3.4.  Modelling  of  combustor  and  other  components 


It  is  assumed  that  the  exhaust  gas  from  the  SOFC  is  com¬ 
pletely  burnt  in  the  combustors,  i.e.,  CO  and  H2  in  the  exhaust 
gas  are  converted  to  CO2  and  H2O,  respectively.  Thus,  the  heat- 
balance  equation  of  the  combustor  is: 


2>;a«; = 


J 


(17) 


The  heat-balance  equation  for  the  heat  exchangers  is  similar  to 
the  above  equation  with  no  species  change  in  the  heat  exchang¬ 
ers. 


4.  Calculations  for  system  model 

A  flow  chart  of  the  modelling  calculations  for  the  SOFC-AC 
system  is  presented  in  Fig.  3.  The  results  of  every  block  in  the 
flow  chart  become  the  inputs  for  the  next  block.  After  setting 
the  process  parameters,  the  calculation  block  belonging  to  the 
SOFC  model  will  be  executed  several  times  until  the  error  of  the 
energy  balance  is  less  than  the  assigned  tolerance  (0. 1  K).  Other 
blocks  in  Fig.  3  will  only  be  executed  once.  In  this  study,  the 
SOFC  is  assumed  to  be  in  thermal  equilibrium  with  the  exhaust 
gas.  Thus,  by  making  an  initial  guess  of  the  cell  temperature, 
the  equilibrium  data  for  the  three  reactions  in  the  SOFC  can  be 
obtained.  All  characteristics  of  the  exhaust  gas  from  the  anode 
and  cathode  can  be  obtained  by  solving  Eqs.  (6)-(9).  Using  this 
data,  the  SOFC  performance  can  be  obtained  in  terms  of  cell 
voltage  and  cell  current.  The  energy-balance  equation  can  be 
used  to  assess  the  validity  of  the  guessed  temperature.  When 
the  cell  temperature  is  confirmed  and  the  characteristics  of  the 
exhaust  gas  from  the  SOFC  are  calculated,  the  subsequent  cal¬ 
culation  blocks  will  be  solved  one  by  one.  The  output  data  of 
the  former  block  is  the  input  of  the  next  block.  The  numerical 
method  for  the  calculation  of  adsorption  cycle  block  has  been 
described  previously  [12].  Visual  FORTRAN  software  was  used 
to  develop  the  simulation  code  for  this  SOFC-AC  model. 


Output  results 

Fig.  3.  Flow  chart  of  calculations  for  SOFC-AC  system. 


formance  of  SOFC  systems.  The  operating  parameters  for  the 
base  case  are  shown  in  Table  1 .  In  this  study,  heat  power  is  not 
considered.  Thus,  the  total  efficiency  of  SOFC-AC  system  is 
given  by: 


+  F( 


cw 


LHV 


(18) 


where  Pe w  is  the  electrical  power  delivered  by  the  SOFC  and 
LHV  is  the  lower  heating  value  of  the  inlet  fuel.  Simulated  results 
for  the  base  case  show  that  the  SOFC  operates  at  a  voltage  of 
0.775  V  and  a  current  density  of  150  mA  cm-2.  The  SOFC-AC 
cogeneration  system  can  achieve  a  total  efficiency  of  more  than 
77%,  i.e.  electrical  power  62%  and  cooling  power  15%.  About 
1800kW  of  electrical  power  and  450  kW  of  cooling  power  can 
be  produced  by  the  system.  The  simulated  stream  properties  for 
the  base  case  are  shown  in  Table  2. 


5.  Results  and  discussion 

The  results  of  the  numerical  model  can  be  used  to  carry  out 
a  detailed  parametric  analysis  which  can  provide  an  insight  into 
the  effects  of  variation  of  the  operating  parameters  on  the  per¬ 


5.1.  Effect  of  inlet  fuel  flow -rate 

The  effect  of  inlet  fuel  flow-rate  on  system  efficiency  is  pre¬ 
sented  in  Fig.  4.  Here  the  ratio  of  the  inlet  fuel  flow-rate  to  the 
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Table  1 


Prescribed  values  of  parameters  for  base  case 


Parameters 

Value 

Fuel  inlet  composition 

CH4:  84%,  N2:  15%,  C02:  1% 

Fuel  inlet  temperature 

300  K 

Air  inlet  temperature 

300  K 

Air  temperature  after  preheat 

930  K 

Cell  operating  pressure 

1 .0  bar 

dc/ac  inverter  efficiency 

95% 

Burner  efficiency 

100% 

Fuel  utilization  factor,  Uf 

0.8 

Circulation  ratio,  ra 

0.2 

Electrolyte  thickness 

0.015  cm 

Interconnect  thickness 

0.010cm 

Cathode  thickness 

0.20  cm 

Anode  thickness 

0.020  cm 

Half  cycle  time  of  adsorbent  cycle,  />, 

1200  s 

Mass  of  the  adsorbent,  ms 

5  kg 

Number  of  adsorption  units,  N 

250 

Fig.  4.  Effect  of  fuel  flow-rate  on  efficiency. 


inlet  air  fuel-flow  rate  is  fixed.  The  data  show  that  both  the  electri¬ 
cal  efficiency  and  total  efficiency  decrease  as  the  fuel  flow-rate 
increases,  while  the  cooling  efficiency  has  a  maximum  value 
for  the  range  of  fuel  flow-rates  investigated.  The  effect  of  fuel 


Fuel  flow  rate,  kmol/h 

Fig.  5.  Effect  of  fuel  flow-rate  on  cell  voltage  and  cell  temperature. 


flow-rate  on  cell  voltage  and  cell  temperature  is  presented  in 
Fig.  5.  For  a  constant  fuel  utilization  factor  and  circulation  ratio, 
a  higher  rate  of  fuel  flow  means  that  more  current  will  be  pro¬ 
duced  and  hence,  the  current  density  is  increased.  The  higher 
the  current  density,  the  higher  is  the  value  of  the  total  overpo¬ 
tential  produced.  The  increase  in  overpotential  will  result  in  a 
reduction  in  the  voltage  of  SOFC  (see  Fig.  5).  Since  current  has 
a  linear  relationship  with  fuel  flow-rate,  the  electrical  efficiency 
will  decrease.  It  is  shown  in  Fig.  5  that  the  cell  temperature 
increases  linearly  with  an  increase  in  the  inlet  flow-rate.  The 
increase  in  overpotential  will  increase  heat  generation,  which 
will  raise  the  cell  temperature.  The  effect  of  inlet  fuel  flow- 
rate  on  cooling  power  is  shown  in  Fig.  6.  The  cooling  power 
increases  asymptotically  to  a  constant  value  with  increase  in  the 
fuel  flow-rate.  When  the  fuel  flow-rate  increases,  the  electri¬ 
cal  efficiency  will  decrease.  This  means  that  more  heat  can  be 
generated  in  the  combustor.  Thus,  the  driven  temperature  of  the 
adsorption  chiller  also  increases  and  this  leads  to  an  increase 
in  cooling  power  for  a  fixed  cycle  time.  When  the  driven  tem¬ 
perature  increases  to  a  certain  value,  no  more  refrigerant  can 
be  drawn  from  the  adsorbent.  Thus,  the  value  of  the  cooling 
power  will  tend  to  a  constant  value  and  the  cooling  efficiency  will 
decrease. 


Table  2 

Stream  properties  for  SOFC-AC  system 


Temperature  (K) 

Gas  flow  rate  (kmol  h  1 ) 

Gas  composition  (%) 

ch4 

h2 

CO 

co2 

h2o 

n2 

o2 

1 

300 

16 

84.0 

— 

— 

1.0 

— 

15.0 

— 

2 

300 

180 

— 

— 

— 

— 

— 

79 

21 

3 

608 

26.7 

50.3 

5.5 

2.9 

10.4 

19.7 

11.2 

— 

4 

930 

180 

— 

— 

— 

— 

— 

79 

21 

5 

1118 

53.6 

— 

13.6 

7.2 

24.5 

49 

5.7 

— 

6 

1118 

157.6 

— 

— 

— 

— 

— 

90.2 

9.8 

7 

1118 

10.7 

— 

13.6 

7.2 

24.5 

49 

5.7 

— 

8 

1481 

196 

— 

— 

— 

6.9 

13.7 

73.8 

5.6 

9 

934 

196 

— 

— 

— 

6.9 

13.7 

73.8 

5.6 

10 

a 

196 

— 

— 

— 

6.9 

13.7 

73.8 

5.6 

a  Temperature  of  stream  10  will  change  with  time. 
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Fig.  6.  Effect  of  fuel  flow-rate  on  cooling  power  produced. 


5.2.  Effect  of  fuel  utilization  factor 

The  effect  of  the  fuel  utilization  factor  Uf,  on  system  perfor¬ 
mance  is  shown  in  Fig.  7,  while  the  effect  on  cell  voltage  and 
cell  temperature  is  given  in  Fig.  8.  When  the  fuel  utilization  fac¬ 
tor  increases,  the  content  of  hydrogen  in  the  anode  is  reduced. 
A  higher  Uf  will  lead  to  a  higher  value  of  current  and  a  higher 
cell  temperature  (see  Fig.  8).  The  cell  temperature  reduces  the 
overpotential  and  the  current  gives  rise  to  an  increased  over¬ 
potential.  The  SOFC  voltage  will  increase  slightly  and  then 
decrease  quickly  with  increase  Uf  due  to  the  combined  effect  of 
cell  temperature  and  current  density.  The  electrical  efficiency, 
which  is  proportional  to  the  product  of  voltage  and  Uf,  has  a 
maximum  value  at  a  Uf  around  0.85.  For  a  Uf  value  larger  than 
0.85,  the  electrical  efficiency  will  decrease  due  to  the  increased 
effects  of  the  voltage  drop.  A  higher  Uf  also  results  in  a  lower 
combustion  temperature,  that  leads  to  a  lower  value  of  the  driven 
temperature  for  the  adsorption  chiller.  Therefore,  the  cooling 
power  will  decrease  with  higher  Uf.  When  Uf  increases  beyond 


Fig.  8.  Effect  of  fuel  utilization  factor  on  cell  voltage  and  cell  temperature. 

0.85,  the  electrical  efficiency  will  decrease,  which  in  turn  leads 
to  a  higher  depleted  fuel  heating  value.  Thus,  more  heat  will 
be  generated  in  the  combustor  and  the  driven  temperature  of 
adsorption  chiller  will  increase.  Therefore,  the  cooling  power 
and  the  cooling  efficiency  will  increase  with  increase  in  fuel 
utilization  factor  when  Uf  >0.85. 

5.3.  Effect  of  circulation  ratio 

The  effect  of  circulation  ratio  on  system  performance  is 
shown  in  Fig.  9.  It  can  be  seen  that  an  increase  in  the  circulation 
ratio  has  a  positive  impact  on  the  system  electrical  efficiency, 
but  adversely  affects  the  cooling  efficiency.  The  total  efficiency 
can  increase  slightly  with  an  increase  in  the  circulation  ratio.  If 
the  circulation  ratio  increases,  the  global  fuel  utilization  factor 
will  also  increase.  Hence,  the  electrical  efficiency  increases  for  a 
higher  value  of  the  circulation  ratio.  A  high  circulation  ratio  also 
leads  to  a  lower  depleted  fuel  heating  value,  and  hence,  a  lower 
driven  temperature  in  the  adsorption  chiller.  The  latter  leads  to 
a  decrease  in  cooling  power  for  a  fixed  cycle  time.  Hence,  for 
fixed  inlet  fuel  flow-rate,  the  cooling  efficiency  will  decrease 
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Fig.  7.  Effect  of  fuel  utilization  factor  on  efficiency. 


Fig.  9.  Effect  of  circulation  ratio  on  efficiency. 


Y.  Liu,  K.C.  Leong  /  Journal  of  Power  Sources  159  (2006)  501-508 


507 


<D 

V— 

CO 

i_ 

0 

Q_ 

E 

0) 

0 


Fig.  10.  Effect  of  circulation  ratio  on  cell  voltage  and  cell  temperature. 


Inlet  air  temperature,  K 

Fig.  12.  Effect  of  inlet  air  temperature  on  cell  voltage  and  cell  temperature. 


with  increase  in  the  circulation  ratio.  The  effect  of  circulation 
ratio  on  cell  voltage  and  cell  temperature  is  given  in  Fig.  10. 
The  trends  of  the  curves  are  very  similar  to  those  in  Fig.  8  as 
an  increase  of  circulation  ratio  leads  to  a  larger  global  utiliza¬ 
tion  factor.  Thus,  effects  on  cell  voltage  and  cell  temperature  for 
these  two  operating  parameters  are  similar. 

5.4.  Effect  of  inlet  air  preheat  temperature 

The  inlet  air  temperature  can  be  increased  by  increasing  the 
heat-exchange  efficiency  between  air  and  the  mixed  gas  from  the 
combustor  in  heat  exchanger  1 .  The  influence  of  inlet  air  temper¬ 
ature  on  system  efficiency  is  given  in  Fig.  1 1 .  An  increase  in  the 
inlet  air  temperature  yields  an  increase  in  the  electrical  efficiency 
and  total  efficiency  which  have  a  negative  effect  on  cooling  effi¬ 
ciency.  A  high  value  of  the  inlet  air  temperature  will  result  in  a 
high  cell  temperature  and  therefore  result  in  an  increase  in  the 
voltage  (see  Fig.  12).  A  high  inlet  air  temperature  also  leads  to 
a  reduction  in  the  temperature  of  the  steam,  which  drives  the 
adsorber.  Thus,  the  cooling  power  will  be  reduced  with  increase 
in  the  inlet  air  temperature. 


5.5.  Effect  of  mass  of  adsorbent 

The  effect  of  the  mass  of  adsorbent  on  cooling  power  is 
presented  in  Fig.  13.  The  mass  of  zeolite  adsorbent  in  the  adsor¬ 
ber  is  proportional  to  the  number  of  adsorption  units.  There¬ 
fore,  an  increased  mass  of  zeolite  will  lead  to  a  reduction  of 
the  inlet  mass  flow  of  the  exhaust  gas  from  heat  exchanger 
1  for  every  adsorption  unit.  It  was  concluded  in  a  previous 
study  [13]  that  the  thermal  performance  of  the  adsorption  cycle 
will  change  very  little  when  the  mass  flow-rate  of  the  heat- 
exchange  fluid  is  larger  than  a  certain  value.  Hence,  when  the 
total  mass  of  adsorbent  increases,  the  cycled  refrigerant  pro¬ 
duced  by  every  adsorption  unit  changes  negligibly.  Thus,  the 
total  cycled  refrigerant  is  increased  for  a  larger  mass  of  adsor¬ 
ber,  which  leads  to  a  higher  value  of  cooling  power.  With  a 
further  increase  in  the  mass  of  the  adsorbent,  the  maximum 
temperatures  of  adsorbent  decreases  rapidly  and  causes  a  sig¬ 
nificant  reduction  in  the  cycled  mass.  Accordingly,  the  total 
cycled  refrigerant  begins  to  decrease.  In  summary,  the  cool¬ 
ing  power  has  a  maximum  value  with  variation  of  the  mass  of 
adsorbent. 


Fig.  13.  Effect  of  mass  of  adsorbent  on  cooling  power  produced. 
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6.  Conclusions 

The  effects  of  fuel  flow-rate,  circulation  ratio,  fuel  utilization 
factor,  inlet  air  temperature  and  mass  of  adsorbent  on  the  perfor¬ 
mance  of  a  IRSOC-absorption  chiller  system  are  investigated. 
Based  on  simulation  results,  the  following  conclusions  can  be 
drawn. 

( 1 )  The  proposed  IRSOFC-adsorption  chiller  cogeneration  sys¬ 
tem  can  achieve  a  total  efficiency  of  more  than  77%  (elec¬ 
trical  power  of  62%  and  cooling  power  of  15%). 

(2)  Both  the  electrical  efficiency  and  total  efficiency  decrease  as 
the  fuel  flow-rate  increases.  In  addition,  the  cooling  power 
increases  asymptotically  to  a  constant  value  with  increase 
in  the  fuel  flow-rate. 

(3)  The  electrical  efficiency  has  a  maximum  value  at  a  value  of 
Uf  of  0.85  within  the  investigated  range. 

(4)  An  increase  in  the  circulation  ratio  has  a  positive  impact 
on  the  system  electrical  efficiency  but  adversely  affects  the 
cooling  power. 

(5)  An  increase  in  the  inlet  air  temperature  yields  an  increase 
in  both  the  electrical  efficiency  and  the  total  efficiency.  The 


cooling  power  exhibits  a  maximum  value  with  variation  of 
the  mass  of  adsorbent  in  the  adsorption  chiller. 

References 

[1]  J.L.  Silveria,  E.M.  Leal,  L.  Ragonha  Jr.,  Energy  26  (2001)  891-904. 

[2]  M.  Burer,  K.  Tanaka,  D.  Favrat,  K.  Yamada,  Energy  28  (2003)  497- 
518. 

[3]  N.  Douss,  F.  Meunier,  Chem.  Eng.  Sci.  44  (1989)  225-235. 

[4]  A.F.  Massardo,  F.  Lubelli,  J.  Eng.  Gas  Turbines  Power  122  (2000)  27- 
35. 

[5]  S.  Campanari,  J.  Power  Sources  92  (2001)  26-34. 

[6]  S.H.  Chan,  H.K.  Ho,  Y.  Tian,  J.  Power  Sources  109  (2002)  111-120. 

[7]  W.  Zhang,  E.  Croiset,  PL.  Douglas,  M.W.  Fowler,  E.  Entchev,  Energy  Conv. 
Manage.  46  (2005)  181-196. 

[8]  U.G.  Bossel,  Final  Report  on  SOFC  Data  Facts  and  Figures,  Swiss  Federal 
Office  of  Energy,  Berne,  CH,  1992. 

[9]  D .J.  Hall,  Transient  modeling  and  simulation  of  a  solid  oxide  fuel  cell,  PhD 
Thesis,  University  of  Pittsburgh,  1997. 

[10]  R.  Bove,  P.  Lunghi,  N.M.  Sammes,  Int.  J.  Hydrogen  Energy  30  (2005) 
181-187. 

[11]  N.F.  Bessette,  Modeling  and  simulation  of  a  solid  oxide  fuel  cell  power 
system,  PhD  Thesis,  Georgia  Institute  of  Technology,  1994. 

[12]  K.C.  Leong,  Y.  Liu,  Appl.  Therm.  Eng.  24  (2004)  2359-2374. 

[13]  Y.  Liu,  K.C.  Leong,  Appl.  Therm.  Eng.  25  (2005)  1403-1418. 


